home *** CD-ROM | disk | FTP | other *** search
/ Loadstar 179 / 179.d81 / binary stars (.txt) < prev    next >
Commodore BASIC  |  2022-08-26  |  2KB  |  74 lines
  1. 5 poke55,.:poke56,56:clr
  2. 10 dv=peek(186):ifdv<8thendv=8
  3. 15 poke53280,0:poke53281,0:print"[147]"
  4. 16 poke53371,0
  5. 25 poke53272,31
  6. 30 ad=49152
  7. 32 sysad:sysad+12
  8. 35 p1=4*atn(1):r1=180/p1
  9. 38 tq=ti
  10. 40 c=2*p1:a1=.0000005
  11. 42 bx$="[153][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180]"
  12. 44 print"[147][153][219][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][220]":sysad+9,4
  13. 46 printbx$""tab(38)bx$
  14. 48 print"[153][255][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][168]"
  15. 50 print""tab(1)"[129]    [194][201][206][193][210][217] [211][212][193][210] [207][210][194][201][212][193][204] [195][193][204][195][213][204][193][212][207][210]   "
  16. 55 printtab(4)"[156][211]emimajor [193]xis [158][219]s: ";:l9%=6:gosub900:a=q9
  17. 60 printtab(4)"[156][197]ccentricity [158](0-.99) [219]s: ";:l9%=6:gosub900:e0=q9
  18. 65 ifint(e0)=0then80
  19. 70 ife0=>1thenprint"[145][145]":goto60
  20. 75 :
  21. 80 printtab(3)"[156][197]poch of [208]eriastron [158][219]s: ";:l9%=6:gosub900:t=q9
  22. 85 printtab(3)"[156][207]rbital [208]eriod [158]([217]ears): ";:l9%=3:gosub900:p=q9
  23. 90 printtab(4)"[156][201]nclination [158][219]s: ";:l9%=6:gosub900:i=q9
  24. 95 printtab(3)"[156][193]rg. of [208]eriastron [158][219]: ";:l9%=6:gosub900:w=q9
  25. 100 printtab(3)"[156][208].[193]. of [193]sc [206]ode [158][219]: ";:l9%=6:gosub900:l=q9
  26. 110 i=i/r1:w=w/r1:l=l/r1
  27. 115 :
  28. 120 printtab(3)"[156][196]ate of [207]bservation [158]([217]rs): ";:l9%=6:gosub900:d=q9
  29. 125 t1=d-t:m=c*t1/p
  30. 130 e=m+e0*sin(m)+e0*e0*sin(2*m)/2
  31. 135 n=e-e0*sin(e):g=m-n
  32. 140 f=g/(1-e0*cos(e)):e=e+f
  33. 145 ifabs(g)>a1then135
  34. 150 u=(1+e0)/(1-e0):h=cos(e/2)
  35. 155 ifh=0thenv=p1:goto165
  36. 160 v=2*atn(sqr(u)*sin(e/2)/h)
  37. 165 k=v+w:y=sin(k)*cos(i)
  38. 170 x=cos(k):q=atn(y/x)
  39. 175 ifx<0thenq=q+p1
  40. 180 ifx<0then190
  41. 185 ify<0thenq=q+c
  42. 190 p2=q+l:ifp2>cthenp2=p2-c
  43. 195 r=a-a*e0*cos(e)
  44. 200 s=r*x/cos(q)
  45. 205 p3=int(p2*r1*10+.5)/10
  46. 210 s3=int(s*100+.5)/100
  47. 212 printtab(1)"[153][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175][175]":sysad+9,6
  48. 214 printtab(7)"[158][211][197][195][207][206][196][193][210][217] [211][212][193][210][211]-"
  49. 215 printtab(5)"[154][208]osition [193]ngle. ";p3;"[219]"
  50. 220 printtab(5)"[154][211]eparation ";s3;"[154][193]rc [211]econds."
  51. 230 gosub3000
  52. 235 goto44
  53. 900 q9$="":poke198,.
  54. 905 geta$
  55. 910 poke646,rnd(1)*15+1:print"*[157]";:ifa$=""then905
  56. 915 ifa$=chr$(13)thenprint" ":q9=val(q9$):return
  57. 920 if(a$=chr$(20)andlen(q9$))thenq9$=left$(q9$,len(q9$)-1):goto950
  58. 925 iflen(q9$)>=l9%thensysad+9,20:goto905
  59. 930 if(a$>="0"anda$<="9")ora$="."then940
  60. 935 goto905
  61. 940 q9$=q9$+a$
  62. 945 print""a$;:sysad+9,5:goto905
  63. 950 print" [157][157] [157]";:goto905
  64. 3000 poke214,20:print:printtab(8)"[159](1[159]) [212]ry another one
  65. 3010 [153][163]8)"open(2open) (NULL)o (NULL)(NULL)right$(NULL)val(NULL)(NULL)val (NULL)enu
  66. 3020 poke198,0
  67. 3030 geta$:ifa$<"1"ora$>"2"then3030
  68. 3040 ifa$="1"thenreturn
  69. 3050 sysad+15
  70. 3060 print"[147][144]load"chr$(34)"b.universe iv"chr$(34)","dv
  71. 3070 print"run28"
  72. 3080 poke631,13:poke632,13:poke198,2:end
  73. 10000 d=peek(186):n$="binary stars":open15,d,15,"s0:"+n$:close15:saven$,d:end
  74.